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ABSTRACT 

Nuclcosynthctic yield predictions for multi-dimensional simulations of thermonuclear 
supernovae generally rely on the tracer particle method to obtain isotopic information 
of the ejected material for a given supernova simulation. We investigate how many 
tracer particles are required to determine converged integrated total nucleosynthetic 
yields. For this purpose, we conduct a resolution study in the number of tracer particles 
for different hydrodynamical explosion models at fixed spatial resolution. We perform 
hydrodynamic simulations on a co-expanding Eulcrian grid in two dimensions assum- 
ing rotational symmetry for both pure deflagration and delayed detonation Type la 
supernova explosions. Within a given explosion model, we vary the number of tracer 
particles to determine the minimum needed for the method to give a robust prediction 
of the integrated yields of the most abundant nuclides. For the first time, we relax 
the usual assumption of constant tracer particle mass and introduce a radially vary- 
ing distribution of tracer particle masses. We find that the nucleosynthetic yields of 
the most abundant species (mass fraction > 10~ 5 ) are reasonably well predicted for 
a tracer number as small as 32 per axis and direction - more or less independent of 
the explosion model. We conclude that the number of tracer particles that were used 
in extant published works appear to have been sufficient as far as integrated yields 
are concerned for the most copiously produced nuclides. Additionally we find that a 
suitably chosen tracer mass distribution can improve convergence for nuclei produced 
in the outer layer of the supernova where the constant tracer mass prescription suffers 
from poor spatial resolution. 

Key words: nuclear reactions, nucleosynthesis, abundances - supernovae: general 



1 INTRODUCTION 

The basic processes of Type la supernovae (SNe la) 
have been proposed a lmost 50 years ago (see e.g. 
iHillebrandt fe Niemeverl 120001 . for a review): A ther- 
monuclear exp l osion in electron-degenerate matter 
|Hovle fc Fowled Il960l ) produces radioactive 56 Ni that 
by its decay delivers en ergy at exact l y the rate observed 
in SN la light curves (iPankevI fl962l; iTruran et all 1 19671 : 
I Colgate fc McKed Il969l : iKuchner et all 119941 ). Despite 



IWooslev et al.l 120091 ). then better agree ment of the models 
with o bservations can be ob t ained ( e.g. iRopke fc Niemeveil 



l2007l : iBravo fc Garcia-Senj I2008J ; iKasen et all 120091 ). If 



such a d eflagration to det onation transition does not occur 
(cf. e.g. iNiemeverl ll999T ). the possibility that off-center 
ignition in a single spot leads to a detonation of a white 
dwarf out of hydrostatic equilibrium in the so -called GCD 
model is yet another proposed mechanism (jPlewa et al.l 



this long history, the questions of the progenitor system 
and the explosion scenario are not completely answered. 
Pure deflagrations in Chandrasekhar-mass white dwarfs, 
currently somewhat disfavoured for their apparent in- 
ability to produce bright explosions, have been one of 
the contending explosion scenarios for a long time (e.; 
iNomoto et al||l984l ; iGamezo et aHl2003l ; IRopke et ai1l2007 
If th e initial deflagra t ion c a n transition i n to a detonation 
(e.g. iKhokhlov etaH Il997l : iRopkd 120071 : IWooslevI 120071 : 



12004 IJordan et al.|[200a : lMeakin et alj|2009l ). although the 
robustness of the initiati on of the detonation in this model is 
also fa r from certain (cf. IRopke et a,l.ll2007l : ISeitenzahl et al.l 
l2009al lbl). However, there are indicat ions from recent 
stellar population synt hesis studies (e.g. iRuiter et al.ll2009l ; 
iMennekens et al.ll2010l ) and X-ray observations of ell iptical 
galaxies and galaxy bulges IjGilfanov fc Bogdanll2010l ). that 
the long favoured single degenerate Chandrasekhar-mass 
progenitor chan nel is unable to accou nt for observational 
rates of SNe la (jCappellaro et al.lll999l ). Sup ernovae result- 
ing from the merger of two white dwarfs (e.g. IPakmor et al.l 
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[2010J) have more favourable statistics and remain a possible 
explosion channel. Last but not least, the double detona- 
tion sub-Chandrasekhar mass model, which has received 
renew ed interest of late (JFink et alj|2007l , l2010l : ISim et all 
l20ld ), is currently again considered a serious alternative. 

The different nucleosynthesis occurring in all these ex- 
plosion scenarios is an important test for the validity of the 
respective models. For one-dimensional simulations, the nu- 
cleosynthesis can be calculated during the simulation via 
the coupling of a nuclear reaction network to the hydro- 
dynamics. Since today's most promising explosion mod- 
els either explicitly break spherical symmetry (e.g. merg- 
ers, double detonations, off-center ignitions) or, in the case 
of centrally ignited spherically symmetric explosions in- 
clude buoyancy driven turbulent combustion, at least two- 
dimensional simulations are required to simulate the essen- 
tial physics of the explosion. Unfortunately, in highly re- 
solved two-dimensional simulations and especially in three- 
dimensional simulations use of an extended, full nuclear re- 
action network during the hydrodynamic evolution is com- 
putationally not feasible. Consequently, the nuclear energy 
release for multi-dimensional simulations is modeled with 
simplified and approximate schemes (e.g . Khokhlovl 199a; 



Reinecke et al.l 2002al; IVladimirova et 



20071 ; iTownslev et~aT1 120071 ; Scitcnza hTet alj 12009 J ). Al 



I e.g. lAti 
aQboM 



Calder et al.l 



though some works on multi-dimensional explosion sim- 
ulations of SN la do not calculate any de t ailed nucle- 
osynthesis at all (e . g. I Reinecke et al. 200 2bl; IPlewa et al.l 



2004 iGamezo et all 120051; iPlewal 1 2007 

Bravo fc Garcia-Senzll2008l ). a detailed isotopic composition 



J. 
)7; 



Jordan et al.1 |200S| : 



of the ejecta is often determined in a post-processing step 
dTrayaglio et al.ll2004 | 2005l; iBrown et al.ll2005l ; iRopke et al.l 
120061 : iFink et aljl20icl : iMaeda et al.ll2010h . While alterna- 
tive methods for obtaining isotopic nucleosy nthetic yields for 
SNe l a are currently und er investigation (|Seitenzahl et al.l 
120081 : iMeakin et aD l2009l 'l . the tracer particle method still 
remains currently the only viable choice for the task. In 
Section [2] we briefly review the tracer particle method and 
introduce how we distribute tracer particles of variable mass 
in our initial models. In Section [3] we give details about our 
hydrodynamic simulations. Section [4] presents the results of 
a convergence study in the number of tracer particles and 
highlights the advantages of using a radially varying dis- 
tribution of tracer masses. We conclude the paper with a 
summary and discussion in Section [5] 



2 TRACER PARTICLE METHOD 

Tracer particles are a Lagrangian component in an Eulc- 
rian grid code. The particles are assigned masses and po- 
sitioned in such a way that a density profile reconstructed 
from their distribution resembles that of the underlying star. 
During the hydrodynamical simulation, they are advected by 
the flow, recording the history of thermodynamic conditions 
along their path. The tracers are nevertheless considered 
massless in the sense that the mass they represent does not 
couple to the hydrodynamic flow via gravity or inertia - they 
are simply passively advected by the flow along streamlines. 
The tracer pa rticle method in multi ple dimensions was 
employed first by iNagataki et al.l <I1997T) in the context o f 
core collapse supernovae (note that lThielemann et al.l (|l986T ) 



already post-processed Lagrangian mass zones of ID hydro- 
dynamical calculations of SNe la). 



2.1 Tracer particle masses and placement 

Historically, N tracer particles are placed in a star of radius 
R and mass M, such that each tracer particle represents the 
same amount of mass m = M/N . However, other choices 
for the distribution of tracer masses are possible and may 
be sometimes preferable. In this paper we consider a tracer 
mass distribution that varies smoothly with radius and ad- 
ditionally fulfills the following criteria: 



for < r < Ri 

for _Ri < r < R 2 

for R 2 < r < R 

where < Ri < R2 < R- 



const, mass per particle 
const, volume per particle (1) 
const, mass per particle, 



This choice of the functional form of the radial distribu- 
tion of the tracer particles is motivated by the following 
considerations: in the inner part of the star (0 < r < R\) 
the burning occurs at high density and tracer particles of 
constant mass are sufficiently spatially dense to resolve the 
abundance gradients. In a transition zone {R\ < r < R2), 
where the density is lower and nuclear burning is incom- 
plete, intermediate mass elements are synthesized. Due to 
the relatively low density, a constant tracer mass approach 
results in a small spatial density of particles there. The con- 
stant volume requirement effectively moves tracers from the 
well sampled inner regions to regions of lower initial density 
where interesting nucleosynthesis occurs. Finally, a constant 
volume approach all the way out to the surface of the star 
would result in a wastefully large amount of tracer particles 
placed into the very outer layers at very low density where 
usually no nucleosynthesis occurs. The transition back to 
constant tracer masses in the very outer layers where hardly 
any nucleosynthesis occurs (R2 < r < R) alleviates this 
problem. 

In the case of equally massive tracer particles, the white 
dwarf is divided into N r shells of equal mass. Then the same 
number of tracers, Ng, is distributed uniformly in mass into 
each of these shells. If the tracers should represent the same 
volume, the shells are chosen to have the same volume in- 
stead. The distribution of the tracers in these shells is done 
uniformly in mass as before. If both approaches are mixed as 
described above, we require in addition that the tracers on 
either side of the interface between the areas of equal-mass 
and equal-volume all have the same mass. Together with the 
constraints |T|), this leads to the following linear equations 
constraining the particle numbers Ni,N2, N3 in each of the 
three regions: 



Mi 

M - Ah 

iVi + iV 2 + N 3 



Vpi 

N 2 

VP2 

N 2 



C2) 



where Mi = M(Ri) is the enclosed mass as a function of 
radius, pt = p(Ri) is the mass density at position Ri and 
V = [MR2 ~ fl?)]/3. 
Note that for a cold WD in hydrostatic equilibrium of known 
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Figure 1. Shown are the mass (dashed lines) and volume (solid 
lines) represented by a tracer particle as a function of its initial 
position in the mass coordinate for N = 256 X 512 tracers. Thin 
(red) lines are for constant tracer particle mass. Thick (black) 
lines are for variable tracer particle mass for M\ = 1.05 Mq and 
M 2 = 1.355 M Q . 



Nuclear burning is included applying a simplified scheme 
|Reinecke et alJIJOOJal ; iFink et alj|20ld ): level sets are used 
to propagate the nuclear burning flames at the correct speed. 
In this thin flame approximation, an immediate energy re- 
lease is performed behind the leve l set representing t he flame 
surface using the new tables from lFink et alj (J2010|) for det- 
onations. A table for deflagrations was calculated similar to 
that work. 

All simulations presented here were performed assuming 
2D rotational symmetry on a grid with cylindrical coordi- 
nates. The initial stellar model was a cold, isothermal (T = 
5 x 10 5 K) white dwarf in hydrostatic equilibrium of mass 
Mwd = 1.401 M Q and electron fraction Y e — 0.49886, which 
corresponds to a central density of p c = 2.9 x 10 9 g cm -3 . For 
nuclear energy generation purposes, the composition was as- 
sumed to be 50% le O and 50% 12 C by mass homogeneously 
throughout the star. The grid resolution was 512 x 512 cells 
for simulations that were restricted to one hemisphere (as- 
suming mirror symmetry across the equator), and 512 x 1024 
cells for simulations that included both hemispheres respec- 
tively. Three different sets of explosion models were investi- 
gated. They are discussed briefly in turn below. 



composition and total mass M, given any one of the quan- 
tities {Ri, Mi, pi} the other two are uniquely determined by 
the equation of state. 
The solution to the set of linear equations ([2j is given by: 

* = Tn*+i N ' (3) 

As = -: N r , 

A + B + 1 ' 

where A = M 1 /(Vp 1 ) and B = (M - M 2 )/(Vp 2 ). 
The tracers are placed into the star according to the rules 
outlined above. An additional small offset is added to the 
coordinates such that each particle has a random initial po- 
sition within its corresponding fluid element. A comparison 
of the constant and variable tracer mass distributions (for 
the same number of total tracer particles) is shown in Fig. [T] 
It is evident from the figure that, for this particular choice 
of Mi = 1.05 Mq and M 2 = 1.355 Mq, in the outer regions 
of the star (M(t)/Mq >1.15) the variable tracer mass ap- 
proach results in a smaller volume represented by each par- 
ticle. This increased spatial resolution at low density comes 
at the expense of slightly larger particle masses in the inner 
regions of the star. 



3 SIMULATIONS 

The code used to simulate the supernova explo s ions is the 
MPA SN la code (see, e.g. jReinecke et al]|l999l . l2002al ). In 
this Eulerian hydrodynamics code the reactive Euler equa- 
tions are solved using a fi nite volume s c heme based on the 
PROMETHEUS code by iFrvxell et all (|l989T ) which is an 
im plementation of the "p i ecewi se parabolic method" (PPM) 
of IColella fc Woodward! () 19841 ). In order to track the ex- 
pan ding WD during explosion, a co-e xpanding unif orm grid 



as in lRopke fc Hillebrandtl (|2005J ), and lRopkd (|2005l ) is used 



3.1 Centrally ignited pure deflagration 

In the pure turbulent deflagration model, the burning was ig- 
nited centrally in a simple spherical shape of radius 150 km 
with a superposed two-period cosine-wave perturbation of 
a mplitude 30 km (the C 3 ignition configuration as described 
bv lReinecke et al.ll 19991 ). This setup assumed mirror symme- 
try across the equator. This explosion produced ~ 0.45 Mq 
of iron group elements (IGE). 



3.2 Centrally ignited delayed detonation 

In the second model, an initial deflagration was ignited cen- 
trally in exactly the same C3 configuration as in the pure 
deflagration runs (see above). This time, the deflagration 
transitioned to a detonation after the flame had entered the 
distributed bur ning regime and t he Karlovitz numbeijj ex- 
ceeded 250 (see lKasen et al.ll2009i ). This setup also assumed 
mirror symmetry across the equator. This explosion pro- 
duced ~0.57M Q of IGE. 



3.3 Multi point ignition delayed detonation 

For the third set of explosion simulations, the deflagration 
was ignited in 100 ignition kernels of radius 6 km randomly 
drawn from a Gaussian radial distribution with a standard 
deviation of 150 km and from a uniform distribution in an- 
gle, resulting in ignition kernels distributed within the in- 
ner 306 km of the WD core (correspondi ng to the ignition 
model DD2DJso_06 of lKasen etaJIl2009l ). The detonation 
was triggered in the same way as the centrally ignited de- 
layed detonation described above. This setup simulated both 



1 The use of the Karlovitz number is not rigorous here as it de- 
pends on the concept of a laminar flame speed which does not 
exist in the distributed burning regime. Nonetheless, it is formally 
used in order to characterize the strength of turbulence. 
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hemispheres independently and represented a full star in ax- 
isymmetry. This explosion produced ~ 1.03 Mq of IGE. 



3.4 Post-processing 

In the post-processing step, isotopic nucleosynthetic yields 
for 384 nuclides are calculated by integration of a nuclear 
reaction network over the recorded temperature and density 
time histories along the paths the tracer particles took in the 
hydrodynamic simulation. Apart from the added feature to 
allow for variable tracer particle masses (see section 12. ip . 
the post-processi ng scheme is identical to the one used in 
iFink et al.l (J2010|), which in turn is based on lTravaglio et al.l 
( 2004T) . Further detailed description of h ow the code solves 
the nuclear network can be found in iThielemann et al.l 



(|l99(t l. [Th"ielemann et"al] (|l996h . and llwamoto et al.l (|l99Sh 



For the re action rate libraries we use th e 2009 release of 
REACLIB (|Rauscher fc T hielemann 20 0lj), and for the wea k 
reaction rates we use ( Langanke fc Martinez-Pinedd 120001 ) . 
For post-processing purposes, it was assumed that the initial 
composition of the unburned material was a homogeneous 
mix consisting of 50% 16 0, 47.5% 12 C, and 2.5% 22 Ne by 
mass (to account for solar metalicity of the ZAMS progeni- 
tor). 



4 CONVERGENCE STUDY 

An important question to ask is how many tracer particles 
are needed to converge on the final integrated mass frac- 
tions after freeze-out. Only a limit ed amount of work ha s 
been done so far in this direction. iTravaglio et al.l (| 20041 ) 
performed in the context of core collapse supernova mod- 
els a ID resolution study on the number of tracer particles 
and conclude that for 2000 zones convergence is reached at 
1000 particles. They stated that the results may not be ap- 
plicable in the multi-dimensional case and that a resolu- 
tion study for 2D is in preparation, which, however, appar- 
ently was not published. iBrown et al.l (|2005r ) post-processed 
a 2D pure deflagration explosion with 10000 (constant mass) 
tracer particles. They found, that when choosing a random 
subset of 5000 tracer particles, the 56 Ni mass changed by 
8% (they did n't state whe t her it increased or decreased). 
Most recently iRopke et al.l (|2006T ). post-processed 3D pure 
deflagration models with 27 3 tracer particles in the produc- 
tions runs. They looked into variations of the 56 Ni mass by 
increasing the number of tracers to 35 3 and concluded that 
differences were on the percent level. To date, a more de- 
tailed study spanning a larger range of tracer particle num- 
ber and including isotopes different from 56 Ni has not been 
done. 

Here we present a more thorough resolution study in 
the number of tracer particles for the three different ex- 
plosion models outlined in Section [3] For the pure defla- 
gration (see Section 13. 1[) and the centrally ignited delayed 
detonation simulations (see Section [3. 2\ we have calculated 
nucleosynthetic yields based on 64, 256, 1024, 4096, 16384 
and 65536 constant mass tracer particles. For the multi- 
point ignition delayed detonation (see Section l3.3p . which 
did not assume mirror symmetry across the equator, we 
have calculated nucleosynthetic yields based on 128, 512, 



2048, 8192, 32768 and 131072 constant mass tracer parti- 
cles. This corresponds in all cases to 8, 16, 32, 64, 128, and 
256 tracer particles per axis and direction. In addition, we 
have calculated yields for the same multi-point ignition de- 
layed detonation model with variable tracer masses using 
Mi = 1.05 M Q and M 2 = 1.355 M Q as introduced in Sec- 
tion [271] For the variable tracer mass case the total numbers 
of tracers were slightly different, i.e. 122, 499, 2024, 8140, 
32664, and 130865. For practical purposes these numbers 
can be considered to be the same resolution. 

In this paper we are less concerned with the exact mag- 
nitude of the mass fraction of a given nuclide, but rather 
with how well the yields of the run that included the most 
tracer particles are reproduced by a run with fewer tracer 
particles and how quickly the various mass fractions con- 
verge with total tracer particle number. 

In Figs. [2] and [3] we show, for sequences of runs with 
increasing tracer particle numbers, percent differences of nu- 
clide mass fractions compared to the mass fractions obtained 
in the run that included the highest number of tracer parti- 
cles, which serve as "reference values". To visually empha- 
size nuclides with larger mass fractions, the symbol radius, 
Si, increases logarithmically from Si =0.1 for nuclides with 
mass fraction Xi ^ 10 -5 to a possible maximum of Si — 30 
in arbitrary units. 

The left column of Fig. [2] shows results for the pure 
deflagration setup, which used tracer particles of constant 
mass (see Section 13 . 1 P . The two cases containing the least 
number of tracer particles, 8 2 and 16 2 , show globally poor 
agreement of the nuclide mass fractions with the reference 
values. The case containing 32 2 particles already shows good 
agreement of most Fe-peak nuclei, but many of the most 
abundant intermediate mass nuclei are under-produced on 
the 5%-level. Agreement with the reference values is only 
slightly improved by going to 64 2 tracer particles. Using 128 2 
tracer particles finally reproduces the reference values of the 
most abundant nuclei globally on the 2%-level, with a few 
exceptions, most notably isotopes of Ne, Mg, and Al. 

The right column of Fig.[2]shows results for the centrally 
ignited delayed detonation setup, which also used tracer par- 
ticles of constant mass (see Section [3]2}. Similar to the pure 
deflagration case, the two low resolution cases show globally 
poor agreement of the nuclide mass fractions with the refer- 
ence values. Already the 32 2 tracer particle case shows rather 
good agreement with the high resolution reference run for 
the majority of the more abundant Fe-peak and intermedi- 
ate mass nuclides. Compared to the pure deflagration, this 
delayed detonation agrees better with its references values 
for intermediate tracer particle numbers (i.e. 32 2 and 64 2 ). 
This can be traced back to the fact that the delayed det- 
onation model synthesizes nuclides such as 28 Si, 32 S, 36 Ar 
or 40 Ca copiously during the detonation phase in a large 
volume where even the constant tracer particle mass imple- 
mentation results in a tracer particle number density high 
enough to adequately sample the morphology. Notable ex- 
ceptions to the better agreement are nuclides that are pre- 
dominantly synthesized in the outer layers of the star during 



2 Small departures from exact powers of 2 were necessary to as- 
sure mass conservation and the constraint that tracer particles in 
the intermediate region all represent the same volume. 
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Figure 2. Final (i = 10 s) nuclide mass fraction differences 



in percent for a sequence of increasing total tracer 



X i (N 2 )-X i (256 2 ) 
Xi(256 2 ) 

particle number compared with the highest resolved case containing 256 2 tracer particles. The underlying hydrodynamical simulation is 
the same central ignition pure deflagration for the left column (see Section 13.111 and the same centrally ignited delayed detonation model 
for the right column (see Section li3.2ll . In all cases the traditional constant tracer particle mass approach was used. The radius, Sj, of the 
markers increases with mass fraction Xj according to Sj = max{0.1, 29.9[log 10 (JQ) + 5]/5 + 0.1} in arbitrary units. 
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Figure 3. Final (i = 10 s) nuclide mass fraction differences 



in percent for a sequence of increasing total tracer 



X i (2xN 2 )-X i (2x256 2 ) 
X;(2x256 2 ) 

particle number compared with the highest resolved case containing 2 X 256 2 tracer particles. The underlying hydrodynamical simulation 
is the same multi-spot ignition delayed detonation model (see Section l3.3l l for all cases. The left column is for tracer particles of constant 
mass, whereas for the right column the variable tracer mass approach was used (see Section |2. H and Fig.fTJ. The radius, Sj, of the markers 
increases with mass fraction X{ according to Si = max{0.1, 29.9[log 10 (Xi) + 5]/5 + 0.1} in arbitrary units. 
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the detonation phase, such as e.g. 20 Ne, 26 Mg or 27 A1 (see 
the discussion of the multi-spot ignition delayed detonation 
model further below). 

Fig. |3] shows results for the multi-spot ignition delayed 
detonation setup (see Section |3~T|) . with constant tracer par- 
ticle mass runs in the left column, and variable tracer par- 
ticle mass runs (see Section I2,l|l in the right column. In 
spite of the differences in explosion model and nucleosyn- 
thetic yields, the sequence of the constant tracer mass case 
is qualitatively very similar to the centrally ignited delayed 
detonation (compare right column of Fig. [2] to the left col- 
umn of Fig. |3J . In all but the lowest resolution case, using 
a variable tracer mass approach results in a global improve- 
ment of the agreement with the reference values and im- 
proved convergence characteristics (compare right and left 
columns of Fig. [3J • This improvement is due to the fact that 
the variable tracer mass approach alleviates the problem of 
low tracer particle density in the outer layers of the star in- 
herent to the constant mass approach. Abundance gradients 
of problematic nuclides such as e.g. 20 Ne (see Fig.|4^,), 26 Mg 
(see Fig. |3p) or 27 A1 (see Fig. OJ.), which are produced in 
an incomplete burn during the detonation phase at low den- 
sity in the outer layers of the star, are much better resolved 
with the variable tracer mass implementation (compare right 
and left panels of Figs. [4^,, |4p) and [5^,. For nuclides which 
qualitatively agree equally well with their reference values 
for both constant and variable tracer mass cases, such as 
e.g. 28 Si, the variable tracer mass implementation results in 
a much improved spatial resolution of the abundance dis- 
tribution morphology in the outer regions at the cost of a 
slightly worse resolution in the inner regions (see Fig. [5p). 



5 SUMMARY AND CONCLUSIONS 

We have performed a resolution study in the number of 
tracer particles required to converge on the integrated nu- 
cleosynthetic yield calculations for three different SN la ex- 
plosion models. We find that for all explosion models, total 
nucleosynthetic yields appear to have converged (with few 
exceptions, see Section 2} for abundant nuclides with mass 
fractions > 10 -5 to better than 2% at 128 tracer particles 
per axis and direction. Agreement at better than the 5% 
level is already achieved for many of the most abundant nu- 
clides at particle numbers as low as 32 per axis and direction. 
Based on these results, we extrapolate that pub- 



lished yields from the literature JTrayaglio et al.ll2 004. 2005; 
iBrown et al1l2005l : iRopke et af1l2006l;lMaeda et alJl201Ch are 

based on sufficiently high numbers of tracer particles per axis 
that an adequate to good prediction of the most abundant 
nuclei is given. Furthermore, we have shown that for isotopes 
whose origin is near the surface of the WD, such as e.g. 20 Ne, 
26 Mg or 27 A1, better convergence can be achieved if the con- 
straint of constant tracer mass is relaxed and a choice for the 
tracer particle masses is made that better spatially resolves 
the outer layers. 

The question of how many particles are required to get 
converged spectra and light curves from subsequent radia- 
tive transfer calculations is related. It seems plausible that a 
convergence of the total integrated yields is a necessary but 
not sufficient constraint. For the radiative transfer not only 
the total mass of a given isotope is important but also its lo- 
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Figure 4. Initial spatial distribution of constant (left, N = 8192) 
and variable (right, N = 8140) mass tracer particles colored by 
the final mass fraction of 20 Ne (top panel) and 26 Mg (bottom 
panel) after freeze-out (t = 10 s). 



cation in velocity space. It seems that the method of variable 
mass tracer particles is also promising in this context. Es- 
pecially simulations where the nucleosynthesis is dominated 
by a detonation should greatly benefit from a variable tracer 
mass distribution, since the mass fractions in the inner re- 
gions vary smoothly radially and therefore a de-refinement 
in the inner part is not a high price to pay for increased 
resolution in the outer part. However, the reconstruction of 
an isotopic abundance for every computational cell of the 
radiative transfer calculation from the tracer particle distri- 
bution is not unique. The resultant spectra and light curves 
depend somewhat on particular choice of reconstruction al- 
gorithm used. These questions will likely be addressed in a 
future study. 
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